THE EXACT RENORMALIZATION GROUP IN ASTROPHYSICS 



JOSE GAITE 
Centra de Astrobiologia, CSIC-INTA, 
Carretera Ajalvir, km 4, 28850 Torrejon de Ardoz, Madrid, Spain 

The coarse-graining operation in hydrodynamics is equivalent to a change of scale 
which can be formalized as a renormalization group transformation. In particular, 
its application to the probability distribution of a self-gravitating fluid yields an 
"exact renormalization group equation" of Fokker-Planck type. Since the time 
evolution of that distribution can also be described by a Fokker-Planck equation, 
we propose a connection between both equations, that is, a connection between 
scale and time evolution. We finally remark on the essentially non-perturbative 
nature of astrophysical problems, which suggests that the exact renormalization 
group is the adequate tool for them. 



1 Dynamical equations and fluctuations 

1.1 Astrophysical hydrodynamics 

The usual hydrodynamic equations are the continuity equation, the Euler (or 
Navier-Stokes) equation describing momentum conservation, and the energy equa- 
tion. To account for the gravitational field, we must add its corresponding equa- 
tions, namely, the usual equations for its divergence and curl. In cosmology, we 
must also account for the global expansion, which, in the Newtonian formulation, 
introduces the scale factor a(t) and the Hubble constant H(t)a 

Then we have five equations altogether to describe the evolution of the fields, 
but we assume that the solution of the energy equation is a polytropic equation of 
state P — ng' 1 . So we have four equations: 

^+3Hq+-V-(qv) = 0, (1) 
at a 

dv 1 1 VP 

— +Hv + -v Vv = g , (2) 

at a a g 

V g = ~4nGa(g~ g b ), V x g = 0; (3) 

Qb is the background (average) density. 

Of course, the two equations for the gravitational field can be reduced to a single 
equation for the gravitational potential. Then we have two scalar equations, that 
is, the continuity and the Poisson equation, and one vector equation, namely the 
Euler velocity equation, for two scalar unknowns (the density and the gravitational 
potential) and one vector unknown (the velocity). H, gb and a are known functions 
of time in a definite cosmological model. For astrophysical situations in which the 
cosmological expansion is irrelevant, we may put a = 1 and Qb = 0, H = 0. 

1.2 Coarse graining 

The hydrodynamic equations are macroscopic equations, which can be assumed 
to follow from microscopic Newtonian mechanics of particles through an averaging 
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process called coarse graining. The average density Ql{t) of an element of volume, 
say a box of side L, centered on the point r is the total mass of the particles inside 
divided by the volume L 3 . One defines similarly the velocity field «£,(r), etc. 

We must regard that the coarse graining operation does not completely eliminate 
the irregular nature of the microscopic variables but leaves a fluctuating component, 
which should vanish as L grows. In general, the use of hydrodynamics implies a 
thermodynamic treatment of the element of volume, to which one can add the usual 
thermodynamic fluctuations. However, since the gravitational force is long ranged, 
the separation of scales on which the theory of thermodynamic fluctuations is based 
is useless here. Therefore, the hydrodynamic equations must be supplemented with 
the fluctuations, as usual in hydrodynamics, but with fluctuations that are of a 
particular nature.^ As one scales L up, the coarse-grained fields change, as well as 
the fluctuations. This is the essence of the renormalization group. 

Besides stochastic fluctuations, which appear as noise terms in the hydrody- 
namic equations, one can also consider a probability distribution of initial fields, 
as usually done in cosmology. The fluctuations are determined by the correlation 
functions of the noise or the initial conditions. 

1.3 Dynamical evolution 

To solve a set of stochastic partial differential equations (SPDE), one must solve 
before the partial differential equations as if the noise were just an external source, 
and then express the correlation functions of the dynamical fields g and v in terms 
of the correlation functions of the noise (or initial conditions). Of course, the 
solution is of stochastic nature and can also be expressed in terms of a probability 
distribution V[g, v,t]: the correlation functions are given by the corresponding 
functional integrals. The SPDE are therefore equivalent to an evolution equation for 
the probability distribution V . We shall show below how to solve the hydrodynamic 
SPDE and the associated equation for V corresponding to a simplified version of 
the total set of equations. 

2 The exact renormalization group equation 

Apart from their time evolution, the coarse-grained variables change (we may also 
say "evolve") with the coarse-graining length L. ^From a field-theory point of view, 
a smoothing of the "bare" fields g and v is necessary to regularize short-distance 
singularities in functional integrals. It is convenient in our context to define a 
general coarse graining procedure with the help of a window function, that is, a 
function that quickly vanishes outside a neighborhood of the origin of size L; typical 
examples are the sharp-cutoff ("top-hat") window and the Gaussian window. Thus, 
we define, for example, a smoothed density field as the convolution 



a A recent discussion of hydrodynamic equations in the presence of gravitation is given by A. 




(4) 



Dommguez. 
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As L runs, the coarse-grained variables change and so does the probability 
distribution V (for simplicity, here we disregard its time dependence, as well as the 
variable v, and write Vl[q\)'- 

d m . , dhxW L ^ , . dlnW? 8 . m . 1 dWf S 2 ^ , . 

= -^ Vl[s] - "ST ^ ( ^ ]) - 2 "dT (5) 

where Wl is the Fourier transform of the window function. This differential equa- 
tion, describing evolution with the scale L, is the exact renormalization group equa- 
tion for Vl- Its proof is exposed in the appendix. A useful remark is that it is a 
sort of Fokker-Planck equation for L-evolution. 

The exact renormalization group equation relies on Wilson non-perturbative 
philosophy of the renormalization group and, in fact, was proposed by Wilson 
himselfu It has been amply used in high-energy and statistical physics, appearing 
in various forms, essentially equivalent. However, while coarse graining is widely 
used in astrophysical hydrodynamics and the effect of changing the coarse-graining 
length has been discussed on several occasions, the exact renormalization group 
equation does not seem to have appeared before in the astrophysical literature.^] 



3 Dynamical renormalization group 

Dynamical scaling arguments, customary in the physics of surface growth]^ lead 
one to the conclusion that, when the initial condition is forgotten, the dynamical 
variables adopt a scaling form such as, for example, Sql ■= Ql — (q) ~ L a / (t / L z ) , 
where lim^^oo f(x) — 1 and lim x ^of(x) ~ x a l z . In words, the fluctuations grow 
with time as a power law and they eventually reach saturation, in which state they 
depend on the coarse-graining length as another power law. An elementary example 
is the (linear) diffusion equation, in which z = 2. 

The problem is to derive the exponents a and z and the scaling function /. 
Here is where the dynamical RG comes into play: similarly to the situation in static 
critical phenomena, in which the possible types of critical behaviour correspond to 
fixed points of the ordinary RG, the possible types of critical dynamics correspond 
to fixed points of the dynamical RG. 

This philosophy is appropriate for the problem on hand, as attested by com- 
puter simulations of the microscopic Newtonian mechanics of particles. Indeed, one 
observes that the fluctuations grow with time. Moreover, if for a fixed time one 
reduces the coarse- graining length, one observes a similar growth of fluctuations. 
It is natural to expect that the asymptotic state corresponds to a fixed point of the 
dynamical RG. Relying on this philosophy, this RG has been employed to determine 
the fixed points and critical exponents of a simplified hydrodynamic self-gravitating 
model.El 

Unfortunately, the usual treatment of the dynamical RG is perturbative. This 
may render it inapplicable for self-gravitating models, as I will argue later. Besides, 
the fixed point structure depends on the properties of the noise. Both problems 

''Nevertheless, it has been noticed that the noint probability distribution Pl{q) '■= T- , l[s( 1 ')} satis-i 
fies a diffusion equation, analogous to Eq. feh, but with only the last term of the right hand sidea 
It was used to find the distribution of collapsed objects with volume L 3 according to their mass. 
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are overcome by the exact renormalization group as we have derived it: it is non- 
perturbative and independent of the noise. However, it may depend on the window 
function. We shall see that this dependence is equivalent to a dependence on the 
noise and enlightens the role of the latter in the dynamics. 



4 Burgers and Kardar-Parisi-Zhang equations 

In cosmology, it is customary to resort to the Zel'dovich (or parallelism) approxima- 
tion, which constrains the velocity to have the same direction that the gravitational 
field has, that is, v(r,t) = F(t) g(r, t).u This approximation greatly simplifies the 
set of hydrodynamic equations, in fact, reducing them to a single equation for the 
(redefined) velocity. Furthermore, in the case that the polytropic index 7 = 2, this 
equation becomes the three-dimensional Burgers equation, giving rise to the adhe- 
sion modell^ If we write the Burgers equation in terms of the velocity (and gravity) 
potential and we add a noise term, we obtain the Kardar-Parisi-Zhang equation, 
originally proposed for the description of surface growth. Its fixed point structure 
corresponding to white noise, as given by the dynamical RG, is welL-knownJj An 
even richer structure is obtained by using power-law correlated noisc.19 

4-1 Fokker-Planck equation for the Kardar-Parisi-Zhang equation. 
The KPZ equation in Fourier space reads 



A 



-TT = ^k(0 q ) + % = fc 2 k + - ^k ■ (k - q)^ k _ q + (6) 

q 

It is just a set of coupled Langevin Eqs. If the "force" F were linear, then the 
equations would decouple and actually reduce to copies of the Langevin Eq. for a 
Brownian oscillator, which is trivial to solve. 

Let us assume Gaussian correlated noise with (rj^t) r/k' {t 1 )) = D(k) <5k,-k' S(t — 
t'). The Fokker-Planck equation for the probability distribution II (0k, t) associated 
to Eqs. © is 

dU d 1 d 2 n 

We can compare it with the exact renormalization group equation, under the as- 
sumption of dynamical scaling and, hence, of equivalence of time and scale evo- 
lution, that is, t <-> InL and II(0k,£) <-> Vl[q\- We see that both equations are 
approximately equivalent, if D(k) «-> L(dPr,/dL) and Fk(0q) is linear. The for- 
mer condition is interesting, for it relates the freedom in the choice of noise with 
the freedom in the choice of window function, relating the fluctuations with the 
coarse-graining operation. The latter condition shows that the nonlinearity tends 
to suppress an exact equivalence. 

c The one-dimensional Burgers equation had been introduced earlier as a toy model for turbulence.0 
Its 3D counterpart represents pressure-less compressible turbulence without vorticity. 
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5 Critical phenomena and fractals 



The traditional theory of critical phenomena, including critical dynamics, applies 
to situations close to thermodynamical equilibrium, in which the fluctuations are 
small. In contrast, here we are dealing with far-from-equilibrium processes that 
lead to very large fluctuations. In fact, the asymptotic state is probably of fractal 
nature. In cosmology, a fractal structure of galaxies is currently discernible on 
small scales. The scale invariance of fractal structures has led to assimilate fractal 
fluctuations with critical fluctuations and, therefore, the scale atpwhich they subside 
(homogeneity scale) with the correlation length. This is wrong.u 

Actually, one can obtain for a particle distribution with two-point correlation 
( r o/ r )~ 2a that the relative fluctuation of the coarse-grained density is given by 

~1 2 \ n 3 /r„\-2a 



«««.)">_ 



p 1 L 3 

where a 3 is the volume per particle, and B a number. We clearly see that, for the 
fractal state, L <C T"o =>■ Sql ^ Q, so that thermodynamics is not applicable. Only 
when L > ro thermodynamics becomes applicable as the distribution becomes 
critical. 

When the fluctuations are very large and actually drive the dynamics, the use 
of perturbation theory is questionable. In this sense, it is also questionable the 
application of the perturbative dynamical RG to the gravitational interaction. One 
should instead employ a non-perturbative method, as the exact renormalization 
group, advocated here^j 
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Appendix: Proof of the exact renormalization group equation 

There are several proofs of the exact ^normalization group equation. Here, we 
adapt the proof provided by T. MorriglH to our context. 

Let us consider the generator of correlations (characteristic function) 

Z[J] := J VqV[q] exp(-J-g) = J Vg exp{-i iLZJ.. G - 1 -^^-V[ eL ]- J-g}. 

(9) 

As we increase L — > L + dL the dynamical variable ql — > Ql+<il, losing de- 
tail; that is, losing high wavenumber modes in the Fourier space. Let us define 
the density QdL '■= QL+dL — Ql, corresponding to the removed high wavenum- 
ber modes. In the Fourier space, the filtering convolution becomes a multiplica^ 
tion: QL(k) = Wi(fe) g(k). Furthermore, we define the filtered power spectrunxl 
P L {k) = G(k) W L {k) 2 . Hence, P L+dL = P L + P dL . 



d ln the more general field of dynamical critical phenomena, other non-perturbative approaches 
have been proposed.!] 
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We decompose the integral for the generator at the scale L + dL into two inte- 
grals: 



e 2 dL 



V[g L + g dL }-J-(g L + g dL )}. (10) 



k L 

Integrating over g d L to get the generator of correlations of gL, 

Z L [ql, J}= f Vg dL cxp{-^- - V[g L + g dL ] - J ■ (g L + g dL )} (11) 

J 2f dL 



1 



exp{-P rfL J - V L [-P dL J + ql]-J- g L } 



(12) 



for a new potential Vl- (To prove it, use the change of variable g d L = g d L+L — gh ) 
To unravel the meaning of Vl, take J(k) — for the modes to be removed (k > 
1/L). Then Zl^l, J] — exp{— Vl[pl] — J • £>l}, showing that it is an effective 
potential. 

/From its definition, 



^Z L [ QL ,J]=- 
Substituting for Zl and operating, 



dV L 
dL 



1 dP L 

2 dL 



fdV L 



\dg L 



1 dP d L 

2 dL 



d 2 V L 

del 



d 



dJ{k) 



QL 



Z L 



1 dP L Vl d 2 

2 dL e dg\ 



(13) 



The first equation is the usual form of the exact renormalization group equation, 
but the second equality allows us to write it as 

1 dP L d 2 



-V L _ 



-Vl 



dL 



2 dL dg\ 



(14) 



Adding the quadratic part to get the full V and substituting for the power spectrum 
in terms of W, we obtain Eq. (§). 
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